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ABSTRACT 

We use a volume- and flux-limited sample of local (0.03 < z < 0.1) radio galaxies 
with optical counterparts to address the question of how long a typical galaxy spends 
in radio-active and quiescent states. The length of the active phase has a strong de- 
pendence on the stellar mass of the host galaxy. Radio sources in the most massive 
hosts are also retriggered more frequently. The time spent in the active phase has the 
same dependence on stellar mass as does the gas cooling rate, suggesting the onset 
of the quiescent phase is due to fuel depletion. We find radio and emission line AGN 
activity to be independent, consistent with these corresponding to different accretion 
states. 
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1 INTRODUCTION 

Recent theoretical and observational evidence suggests that 
the growth of galaxies and supermassive black holes at their 
centres are closely related phenomena (Magorrian et al. 
1998; Gebhardt et al. 2000; Haring & Rix 2004). As accre- 
tion onto these black holes is believed to power Active Galac- 
tic Nuclei (AGN) jets, such AGN activity is thus intricately 
linked to the process of galaxy formation and evolution. The 
notion of AGN feedback through radio sources, where the 
AGN jets heat up and expel the surrounding gas, has re- 
ceived particular attention in recent times for a number of 
reasons. Jet heating of the intracluster gas (ICM) decreases 
the rate of accretion onto the central black hole, until it is 
shut off completely. Once the gas has had sufficient time to 
cool, the accretion can restart. Although one can envisage 
a stable configuration in which the accretion rate and gas 
heating are balanced, outward transport of central gas by the 
jet (e.g. Basson & Alexander 2003) implies that AGN heat- 
ing must be sporadic in order to keep some gas available for 
fuelling (|Kawata fc Gibson 2005 | l . Even if all the fuel is car- 
ried away from the central regions, the AGN jet can restart 
once this gas is replenished, for example via an interaction 
with another galaxy (|Bahcall et al. 1997p . Such an interac- 
tion will eventually lead to an increase in the accretion rate 
and hence a re-launching of the radio jets. This inherent in- 
termittency of the feedback process is confirmed by observa- 
tions of rising radio bubbles (|Churazov et al. 2001[) and mul- 
tiple radio lobe pairs in the same source (Venturi et al. 2004; 
Giovannini et al. 1998, 1999). The intermittency also pro- 
vides a natural way of keeping the black hole and spheroid 
growth in step. Furthermore, AGN heating has also been 



invoked to explain the lack of star formation in the most 
massive galaxies (the so-called cosmic downsizing; e.g. Cro- 
ton et al. 2006), and suppression of cooling flows in the cores 
of massive clusters (|Fabian et al. 2003] IForman et al. 2005P . 

The questions of how long an average radio source 
spends in an active state, and the length of time be- 
tween outbursts, are thus crucial to quantifying the impor- 
tance of AGN feedback. A number of methods have been 
employed in attempts to quantify these. The oldest and 
best-known of these invokes spectral ageing and lobe ex- 
pansion speed arguments ( [Alexander fc Leahy 1987| , yield- 
ing on-times of a few times 10' years. Simulations sug- 
gest the jets will be disrupted after times of the order of 
10** years (Tucker & David 1997; Omma & Binney 2004); 
while energy injection rates required to quench cooling 
flows point to timescales of the order of few times lO''- 
few times 10* years pulsen et al. 2005|[O^en & EilekT998] 
IMcNamara et al. 2005|) . Observations of revived radio relic 
sources ( |Ensslin fc Gopal-Krishna 2001[ ) constrain the time 
between outbursts to around 10^ years. The local radio loud 
fraction of a few percent (e.g. Best et al. 2005b) also sug- 
gests the radio sources spend an order of magnitude longer 
in their quiescent state than in the active state. 

Radio source evolution depends strongly on both 
the enviroment of the radio source, and jet char- 
acteristics (IKaiser fc Alexander 199"7l IKaiser et al. 19971 
[Alexan der '2000] IKaiser fc Cotter 2002)l . Recently, Best 
et al. ()2005b|) have found a strong dependence of the lo- 
cal (0.03 < z < 0.1) radio loud AGN fraction on the stellar 
mass of the host galaxy. In other words, it appears that 
the radio jets get retriggered more frequently in massive 
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hosts. Best et al. (2005a, b) also investigated the relation- 
ship between radio and emission line AGN activity by split- 
ting their radio-optical sample into AGNs and star-forming 
galaxies according to the diagnostics appropriate for each 
case, and found that these are different phenomena. This 
finding is confirmed by studies of radio and emission line 
AGN activity in the brightest group and cluster galaxies 
(|Best et al. 2007) . In this work we investigate these intrigu- 
ing results in more detail by employing a flux- and volume- 
limited sample of radio sources with host properties. 

Much like Best et al. (|2005ap . we combine optical infor- 
mation available from the Sloan Digital Sky Survey (SDSS) 
Data Release 2 (|York et al. 2000"|l with two radio surveys, 
NRAO VLA Sky Survey (NVSS; Condon et al. 1998) and 
Faint Images of the Radio Sky at Twenty Cm (FIRST; 
Becker et al. 1995) to construct a local (0.03 < z < Q.l) 
radio-optical sample. Unlike previous works, which are only 
flux-limited, however, our sample is both volume- and flux- 
limited at both the radio and optical (r-band) wavelengths. 
We also use a combination of the FIRST and NVSS 1.4 GHz 
fluxes, ensuring sensitivity to both compact and diffuse 
structures. 

As the main focus of this work is on radio AGN activ- 
ity, we split our sample into AGNs and star-forming galaxies 
according to the radio diagnostic, and construct the bivari- 
ate radio-optical luminosity function by further dividing the 
AGN subsample in stellar mass. Radio source models are 
then used to fit individual source sizes and luminosities (as 
derived from the FIRST and NVSS catalogues) and derive 
the typical length of the active phase as a function of stellar 
mass. Observed bivariate luminosity functions then set the 
time a radio source spends in a quiescent state. 

The paper is structured as follows. In Section [2] we de- 
scribe our volume- and flux-limited sample. Separation into 
AGNs and star-forming galaxies, as well as the resultant 
radio and bivariate luminosity functions, are discussed in 
Section [31 Section 3] outlines our radio source model, and 
the results are presented in Section [S] We conclude with a 
discussion of our findings in the context of AGN fuelling 
mechanisms in Section [7] 

Throughout the paper we assume the cosmological pa- 
rameters Q,u — 0.3, Q,A = 0.7 and Hq — 70 kms~^ Mpc~^. 



2 SAMPLE 

The sample used in this work was obtained by cross- 
correlating the optical Sloan Digital Sky Survey (SDSS) 
Data Release 2 (|York et al. 2000"1l with two Very Large Ar- 
ray (VLA) radio surveys at 1.4 GHz, the NRAO VLA 
Sky Survey (NVSS; Condon et al. 1998), and Faint Images 
of the Radio Sky at Twenty Centimetres (FIRST; Becker 
et al. 1995). The SDSS catalogue used here derives from 
the volume-limited (0.005 < z < 0.1) sample described by 
Nikolic et al. (|2004p . The sample covers a spectroscopic area 
of 2627 deg^ on the sky and is optically complete, con- 
taining 44,630 galaxies with Petrosian r-band magnitude 
Afr < —20.45. For a complete description of the sample and 
various derived optical parameters the reader is referred to 
Nikolic et al. ((2004)1 . 



2.1 Catalogue construction 

The two radio surveys used in constructing our radio-optical 
sample complement each other, in the sense that NVSS is 
more sensitive to diffuse sources, but provides less accurate 
position estimates (resolution of ~ 45 arcsec compared to 
5 arcsec for FIRST). The process of pairing the radio and 
optical catalogues by comparing celestial positions of indi- 
vidual objects was performed using the Catalogue Utilities 
for Reporting, Selecting and Arithmetic (CURSA) software 
package produced by the Starlink Project. Each of the ra- 
dio surveys was independently paired with SDSS, and the 
resultant catalogues then paired together to eliminate any 
duplicates. 

Crucial to the pairing process is the separation distance 
within which two identifications are assumed to correspond 
to the same object. Two different methods were employed 
to address this issue. In the first method, we followed Ivezic 
et al. (|2002[) and performed the pairing process at various 
separations, ranging from 1-40 arcsec for FIRST/SDSS, and 
1-60 arcsec for NVSS/SDSS. The number of matches was 
plotted as a function of separation distance, and the point 
at which the number of matches stays approximately un- 
changed for any further increase in separation distance was 
taken as the best pairing distance estimate. This yielded 
a separation distance of 3 arcsec for FIRST/SDSS, and 12- 
15 arcsec for NVSS/SDSS. The second method involved gen- 
erating a mock catalogue by offsetting all sources in the 
SDSS catalogue by 4 arcmin in right ascention. Both the 
mock and the actual SDSS catalogues were then paired with 
FIRST and NVSS for a range of separation distances. At the 
largest separation distance, we expect all real matches to be 
captured, plus an unknown number of false matches (due 
to the allowed separation distance being too large). An ap- 
proximate number of real matches was obtained by subtract- 
ing the number of total matches for the real catalogue from 
the corresponding value for the mock catalogue. The sum- 
mary of "true" and "false" matches thus obtained is given 
as a function of separation distance in Table [1] The opti- 
mum separation distance yields the highest possible com- 
pleteness rate (given by the fraction of real matches out 
of all matches), and lowest contamination (fraction of false 
matches). We thus chose separation distances of 5 arcsec 
for the FIRST/SDSS pairing process (giving 2448 matches), 
and 15 arcsec for NVSS/SDSS (giving 2393 matches). 

The two catalogues thus obtained were paired to- 
gether with a separation distance of 0.1 arcsec, the level of 
SDSS coordinate accuracy. There were 1427 sources com- 
mon to both catalogues; 1021 FIRST/SDSS objects did 
not have NVSS/SDSS counterparts, while 966 NVSS/SDSS 
sources were without corresponding identifications in the 
FIRST/SDSS catalogue. 

There are two reasons for this. FIRST covers less of 
the sky than NVSS, meaning some NVSS/SDSS sources will 
not have a corresponding FIRST/SDSS match. The rest of 
the mismatches come about due to the differences in sen- 
sitivity, positional accuracy and resolution of the two radio 
surveys. For example, FIRST has higher resolution, which 
means nearby sources confused in NVSS are often resolved in 
FIRST. Conversely, diffuse radio structures with no bright 
central core are not picked up by FIRST, but appear in 
NVSS. All the mismatches in areas of the sky common to 
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FIRST/SDSS 

Separation Real matches Completeness False matches Contamination 



1.5 arcsec 2243 87.7% 0% 

3 arcsec 2391 93.4% 2 0% 

5 arcsec 2448 95.7% 13 0.5% 

10 arcsec 2553 99.8% 77 3.0% 

20 arcsec 2814 100% 289 11.3% 

30 arcsec 3221 100% 651 25.4% 

40 arcsec 3745 100% 1186 46.3% 



NVSS/SDSS 



Separation 


Real matches 


Completeness 


False matches 


Contamination 


5 arcsec 


1205 


43.5% 


15 


1.2% 


10 arcsec 


2001 


72.2% 


62 


3.1% 


12 arcsec 


2200 


79.4% 


81 


3.7% 


15 arcsec 


2393 


86.3% 


119 


5.0% 


20 arcsec 


2656 


95.8% 


202 


7.6% 


40 arcsec 


3600 


100% 


880 


24.4% 


60 arcsec 


4799 


100% 


2027 


42.2% 



Table 1. FIRST/SDSS and NVSS/SDSS catalogue matching summary. For FIRST/SDSS we adopt 3 arcsec as the separation distance 
with the best compromise between maximum completeness and minimum contamination. For NVSS/SDSS the corresponding value is 
15 arcsec. 



FIRST and NVSS were therefore examined visually, and 
eventually all of these were included in the catalogue, yield- 
ing 3110 sources of NVSS/FIRST/SDSS detections in the 
volume 0.005 <z< 0.1. 

Appropriate radio source fluxes were determined for 
each source using the criteria of Table [2] We note that, in 
contrast to Best et al. (|2005ap . we do not always use the 
NVSS flux for sources which are identified in both radio sur- 
veys. This ensures exclusion of contributions from unrelated 
FIRST sources that appear as a single confused source in the 
NVSS envelope. The catalogue was truncated to exclude all 
sources below the NVSS sensitivity limit of 3.4 mjy (corre- 
sponding to 99% completeness in NVSS, the less sensitive of 
the two radio surveys; Condon et al. 1998). A lower redshift 
limit oi z = 0.03 was imposed in order to ensure accurate 
positional matching of the catalogues. This yielded the final 
volume- and flux-limited sample of 1191 radio sources with 
optical identifications. 



from the centre of the optical galaxy than the NVSS/SDSS 
separation distance of 15 arsec, they are missed by the pair- 
ing process. Hence the presented catalogue does not include 
diffuse (greater than a few tens kpc in the adopted cosmol- 
ogy) radio relics. Best et al. (2005a) constructed a similar 
SDSS/NVSS/FIRST catalogue (albeit using slightly differ- 
ent fiux criteria and without the r-band magnitude cutoff), 
but adopted a more sophisticated pairing process. These au- 
thors used a search radius of 3 arcmin (corresponding to 
a few hundred kpc in the adopted cosmology) to identify 
multi-component NVSS sources. They found 6% of their 
sample to consist of such objects. Since these include sources 
with radio cores that are picked up by the pairing pro- 
cess employed here, missing diffuse sources without compact 
cores is not expected to alter the conclusions of this work sig- 
nificantly. Their relatively small number is consistent with 
the notion that radio relics fade away quickly once energy 
supply to the hotspots ceases (|Giovannini et al. 1999p . 



2.2 Selection effects 

The use of two radio surveys with complementary features 
ensures the sensitivity of the sample to both compact and ex- 
tended structure. FIRST will resolve confused NVSS sources 
into multiple compact components. NVSS, on the other 
hand, picks up diffuse structures that FIRST misses. There 
are, however, some types of sources that the approach out- 
lined above is not sensitive to. One such example includes 
the very extended NVSS sources that would only be identi- 
fied as cores by FIRST/SDSS. This problem is circumvented 
by visual inspection of all 1191 radio sources in the final cat- 
alogue. 

More troublesome are sources with no bright central 
component in either of the two radio surveys. If the extended 
components (such as lobes) of these objects are further away 



3 THE LOCAL RADIO LUMINOSITY 
FUNCTION 

3.1 Observed luminosity function 

We construct the radio luminosity function (RLF) for our 
radio flux- and volume-limited sample. In our cosmology, 
sources with Liacui < 8 X 10^^ WHz-i are too weak to be 
detected throughout the observed volume, and are only seen 
because of their low redshifts. For these sources we apply the 
usual V/Vlnax correction (|Condon 1989|) . The resultant RLF 
is given in Tableland also plotted in Figure ISTT] along with 
relevant previous works. 

The sample is split into star-forming galaxies and 
AGNs. Traditionally, this is done using the emission line 
properties of the objects as given by their location in the 
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Use FIRST Flux 




Use NVSS Flux 


Fl 


source is compact in FIRST 
(possibly extended in NVSS) 
and there are other compact 
sources in the NVSS envelope 
or a source with no obvious 
connection to the central 
FIRST/SDSS source of interest 


Nl 


source is compact in FIRST 
(possibly extended in NVSS) 
with no other obvious unrelated 
sources in the NVSS envelope 


F2 


no catalogued source in NVSS 
because either the NVSS field 
is empty or the fiux detected by 
NVSS is too low to be ofliicially 
identified as a source 


N2 


source is extended in NVSS 
or unresolved and will need 
integrating 


F3 


source is confused in NVSS 
but potentially extended in 
b IRS i and will need mtegratmg 


N3 


source is extended in NVSS 
within the FIRST envelope 
with no other unrelated sources 
in either field and NVSS flux is 
complete 


F4 


calculated integrated fiux is 
higher in FIRST and/or assumed 
to be more accurate 


N4 


source not detected in FIRST 
because emission was too 
diffuse and/or extended 



Table 2. Flux criteria. 



All radio sources radio loud AGNs Star-forming galaxies 

21.90 
22.30 
22.70 
23.10 
23.50 
23.90 
24.30 
24.70 
25.10 
25.50 



Table 3. Local radio luminosity function at 1.4 GHz for all sources, and split up into AGNs and star-forming galaxies using the 4000A 
break classification. Number density $ is in units of Mpc~^ 



[OIII] 5007/H/3-[NlI] 6583/Ha plane l|Baldwin et al. 19811 
ISadler et al. 20021 [Kauffmann et al. 2003bp . However, Best 
et al. (|2005ap found that the radio and emission line AGN 
activity are independent phenomena, and argued that the 
4000 A- break, Dn(4000), represents the best way of separat- 
ing the radio AGNs from their star-forming counterparts. 
This break is a consequence of accumulation of a large num- 
ber of spectral lines in a narrow region in wavelength, and 
departs significantly from unity for older, metal rich galax- 
ies (jKauffmann et al. 2003bjl . Using derived star-formation 
rates. Best et al. (2005a; their Fig. 9) calculated the con- 
tribution of stellar emission to the 1.4 GHz luminosity and 
compared this with measured values to derive a demarcation 
in the Dn(4000)- ^^f^^^ plane separating the objects into 
radio AGNs and star-forming galaxies. Since in this work 
we are primarily concerned with the radio AGN activity, we 
adopt this demarcation. Stellar mass Af* for each object was 
evaluated by using the Petrosian z-band magnitude, and as- 
suming that all stars are of solar type. Corrections for dust 
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obscuration were made by comparing z and K-band fluxes. 
Nikolic ()2005p found z — 7f = 2.5 for our sample, and thus 
the stellar masses derived from z-band luminosities are un- 
derestimated by a factor of 3.2. Proceeding in this way, we 
find 653 radio loud AGNs and 538 star-forming galaxies. The 
resultant RLFs are given in Tableland plotted in Figure (2] 

The total luminosity function is in good agreement 
with previous studies at high luminosities, Li.4ghz > 
lO-^"^ WHz~^. At luminosities below this value, however, 
our luminosity function is systematically lower than those 
of other authors. It is due to the fact that we are using 
a volume-limited sample, and have introduced a cutoff of 
—20.45 in r-band magnitude. The radio-optical correlation 
at low (Li,4GHz < 10^'^ WHz^^) luminosities means that in 
discarding objects with lower r-band magnitudes we are also 
excluding some radio sources. In Figure |3] we illustrate this 
effect by binning our sample in radio luminosity, and plot- 
ting distributions of r-band magnitudes in each bin. There 
are relatively few objects with Mr > —21.8 at higher radio 
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lo'' 10" lo'' lo'"* 10"-' lo'' 10"' 10^' 10" lo" lo" 10"' 10^' 10' 
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(a) AGNs (b) Star-forming galaxies 

Figure 2. Local radio luminosity function at 1.4 GHz for (a) AGNs, and (b) star-forming galaxies as classified by a demarcation in the 
Dn (4000)- -^^J^^Sita. plane l IBest et al. 2005a1l . There is good agreement between our work and other studies for Li.4ghz > 10^^ WHz~^. 
However, the luminosity functions diverge at luminosities below this value. This is especially evident in the case of the star-forming 
galaxies, and is due to our cutoff in r-band optical magnitude and the Mr-I/i,4GHz correlation (see Figurc[3ll. 



this work 
Condon (1989) 




10-' 10-^ 10-' 10-^ 10^-^ 10-" 10-' 

L,4/ WHz"' 

Figure 1. Local radio luminosity function at 1.4 GHz for all 
sources. It is in good agreement with other studies (Sadler et al. 
2002 using 2dFGRS, Machalski & Godlowski 2000 using the 
LCRS, and Condon 1989 using a VLA sample). We also find simi- 
lar good agreement at high (Z/i.4ghz >~ 10^"^ WHz~^) luminosi- 
ties with Best et al.'s I l2005al l RLFs for a flux-Umitcd NVSS/SDSS 
sample. To compare these directly we averaged our luminosity 
functions over one dex in Li.4ghz, instead of the bin width of 0.3 
dex. At Li.4GHz <'~ 10^^ WHz~^ wc significantly underestimate 
the source counts when compared with Best et al. due to our 
sample being volume-limited. 



luminosities; however, when I/i.4ghz < 10^^ WHz-i there 
is a significant contribution from galaxies with Mr > —21.8, 
suggesting that there is also an appreciable number of low 
luminosity radio sources with Mr > —20.45 that we have 
discarded due to our optical completeness limit. Not sur- 
prisingly, the difference between our luminosity function and 
those for flux-limited samples becomes apparent at around 

It is worth noting that the 4000A-break is not the only 
way to separate star-forming galaxies from radio AGNs in 
our sample. Another possible method is to use the concen- 
tration index, defined as the ratio of Petrosian 50-to-90% 
light radii, C — Rbo/Rqo- Low values of C (<~ 0.33) corre- 



22.0 < log Lj 4 / W Hz" < 22.5 ^-m — i 
22.5 < log Lj 4 / W Hz , < 23.0 ^-»-^ 
23.0<logLi4/WHz ' <23.5 ^ * 




Figure 3. Cuts through the bivariate radio-optical luminosity 
function. At high (Li.4ghz > 10^'' WHz"'^) radio luminosities 
only very few objects have low enough r-band absolute mag- 
nitudes to be affected by our optical completeness limit. How- 
ever, objects with optical magnitudes at or below our complete- 
ness cutoff of Mr = —20.45 become important for Z/i.4ghz < 
10^'' WHz^^. These missing sources account for the differences 
between our fiux- and volume-limited luminosity function, and 
those derived from samples that arc only flux-limited. 

spond to early-type galaxies, which are more likely to har- 
bour powerful radio sources, while late-type (star-forming) 
galaxies are those with C > 0.375. We have carried out the 
subsequent analysis using both the Dn(4000) and concen- 
tration index diagnostics to split up our sample, and found 
very similar results. Therefore, only the 4000A-break demar- 
cation findings are presented here. 

3.2 Bivariate luminosity function 

We construct a bivariate luminosity function by binning all 
objects in stellar mass and plotting the radio luminosity 
function for each bin. Fractions plotted represent the num- 
ber of AGNs in a certain stellar mass range brighter than 
a given radio luminosity, divided by the total number of 
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L14/WHZ"' 

Figure 4. Bivariate luminosity function, with objects binned by 
stellar mass as calculated from the z-band magnitude (see text). 
Cumulative fraction of AGNs within a given Mi, range brighter 
than the specified radio luminosity {fraction (L) = — ^ jv '^^^^^ 
is shown. 

objects in that stellar mass bin. Again, corrections using 
the y/V„.ax method (|Condon 1989|l were applied. As r and 
z spectroscopic bands largely track the same population of 
old bulge stars, we expect completeness in r-band to imply 
that our sample is also complete in z-band, and hence in 
stellar mass. The bivariate luminosity function is shown in 
Figure |4l 

The most striking feature of the bivariate luminosity 
function is the strong mass dependence of the radio loud 
fraction (equal to the cumulative fraction in the lowest radio 
luminosity bin). This result was also found by Best et al. 
(|2005b|) . with the radio loud fraction Jrl oc M*^'*. A vital 
question is whether this implies that in the more massive 
hosts the radio sources are on for longer, or whether they 
are simply triggered more frequently (with the duration of 
a typical active phase being the same for all masses). We 
address this question in subsequent sections by employing 
detailed radio source modeling. 



4 RADIO SOURCE MODEL 
4.1 Source evolution 

We model radio loud AGNs as an evolutionary phase in the 
lifetime of every galaxy. The source is radio loud when the 
synchrotron jet is "on" ; and once it switches off or the source 
luminosity falls below the detection threshold, the source be- 
comes radio quiet. We assume the black hole mass is already 
in place when the jet switches on, consistent with results of 
various semianalytic models (e.g. Kauffmann & Haehnelt 
2000; Bower et al. 2006; Croton et al. 2006) and observed 
black hole accretion rates in quasars ( [Hopkins et al. 2006| 
lYu fc Tremaine 200311 . and hence jet injection is represented 
by a top hat function, with the durations of on and off phases 
denoted as to„ and toff- 

To account for the initial rise in the source radio power, 
we assume the source initially evolves in a flat atmosphere 
(such as a galaxy core), followed by two power-law profiles 
of the form p(r) = Pcoro (7^—) corresponding to expan- 



^trans 




Figure 5. Radio source and environmental parameters. The col- 
limated jet expands into the ICM described by a core and double 
power-law, and eventually terminates in a hotspot, inflating a 
cocoon of synchrotron-emitting radio plasma. Because the jet is 
supersonic, a bow shock forms ahead of the cocoon, and a contact 
discontinuity separates the entrained gas from the radio plasma 
(see Kaiser & Alexander 1997). 



sion within a galaxy, followed by a steeper cluster atmo- 
sphere (see Figure [S| . We adopt the models of Alexander 
(|2000|) for the initial evolution within the core, and Kaiser 
& Alexander (|1997p and Kaiser et al. (|1997|l for evolution 
in a power-law profile. As the radio source ages, it will suf- 
fer adiabatic, synchrotron and inverse Compton losses. Once 
energy supply from the jet ceases, we assume the radio lumi- 
nosity quickly drops to a value below our detection thresh- 
old. In practice, the cocoon will enter a "coasting" phase 
(|Kaiser fc Cotter 2002^ : however the radio luminosity de- 
clines rapidly in this phase, rendering our approach suffi- 
ciently accurate. 

These models describe the evolution of powerful (FR- 
II) radio sources. Although our sample consists of both these 
objects and the less powerful FR-Is, our approach is justified 
in the following sections by the relative paucity of resolved 
sources with FR-I morphologies. 

For an active jet, source size and radio luminosity are 
functions of core density Pcorc, core radius rcorc, transition 
radius rtra„a between the galaxy and cluster power laws, den- 
sity exponents /Sgaiaxy and /^duster, jet opening angle 6 (re- 
lated to the axial ratio Rt of the source) and jet power 
Qjct- In adopting values for these parameters we are guided 
by observations of nearby X-ray luminous elliptical galaxies 
and relaxed clusters. In the inner regions ([Alien et al. 2006|) . 
power law exponents /3 ~ 0.8-1.2 and core radius rcopo ~ 
1 kpc provide a good fit to r ~ 10 kpc. Density profiles 
of nearby relajced galaxy clusters (|Vikhlinin et al. 2006[l are 
well-fitted by a double power-law, with inner regions hav- 
ing /Sgaiaxy ~ 0.8-1.1, and /Jciuator ~ 1.8-2.6. Trausitiou be- 
tween the two power laws occurs at rtranu ~ 50-200 kpc. 
As Figure [6] illustrates, radio source tracks in the power- 
size (P-D) plane are rather flat and hence not very sensi- 
tive to the exact value of /Jgaiaxy for 0.8 < /3gaiaxy < 1.2 
until the late phases (1— a few times 10* years) of their life- 
times, when the sources begin to suffer appreciably from in- 
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Figure 6. Example radio luminosity— size tracks for sources 
evolving in different atmospheres. In all cases the adopted pa- 
rameters are Pcoro = 3.7 x 10~^^ kgm~^ (corresponding to 
ricore = 0.2 cm-3), rcora = 1 kpc, /JciuBtor = 1.9, Rt = 2 aud 
Qj3t = 3 X 10''^ W. Evolution for f„„ = 5 X 10* years is shown. 
Varied parameters are the inner density exponent, /3gaiaxy, and 
transition radius rtrana. Long-dashed curve is for /Sgaiaxy = 0.8, 
''trans = 50 kpc; dotted curve is for /3gaiaxy = 1-0, ^trana = 50 kpc; 
short-dashed curve is for /3gaiaxy = 1-2, rt,ana = 50 kpc; and solid 
curve for /Jgaiaxy = 1.0, rtrana = 20 kpc. Also shown are the time 
markers of 10^, lO'', 3 X lO'^, 10* and 3 X 10* years. The source 
lumininosity does not evolve substantially until it begins to suf- 
fer significant inverse Compton losses, which happens around the 
same time as the source enters the steeper part of the atmosphere. 

verse Compton losses (|Kaiser et al. 1997)) and/or enter the 
steeper part of the atmosphere. At that point, the losses 
will typically be the dominant factor in determining the ra- 
dio power of the source. This allows us to fix /3gaiaxy = 1-0, 
/^cluster = 1.9 and rtrana = 50 kpc in our models. Conser- 
vation of mass then sets the electron density within the 
core to be ricoro ~ 0.1— fewtimes 0.1 cm"'', and we adopt 
ricorc = 0.2 cm"^ from the Allen et al. (|2006|l sample. Fol- 
lowing observations of Cygnus A ( |Begelman" fc Ciofli 1989| 
we also set Rt ~ 2.0, corresponding to a jet opening angle 
of 31° (|Kaiser fc Alexander 1997)) . 



4.2 Dependence of predicted RLF on model 
parameters 

Given durations of jet on and off timescales, we can now 
predict the radio luminosity function (RLF) for a popu- 
lation of such sources observed at random stages in their 
evolution. Figure [7] explores the dependence of the RLF on 
input parameters. We adopt the same base quantities as 
above, in addition taking the scatter in jet power (as ex- 
pected from the observed Mbh — Mi, {Haring & Rix 2004 
[Magorrian et al. 1998) ) and Qjct - Mbh (|Allen et al. 2006 
relations) to be 0.8 dex. The jet is active for ton = 2 x 
10* years, and off for time tatt = 2to„. 

To explore the sensitivity of our predicted cumulative 
RLFs on model parameters, each of these parameters is var- 
ied in turn (with the others fixed). The relative unimpor- 
tance of the exact density profile of the atmosphere into 
which the source is expanding, as discussed above, is the 
reason very similar RLFs are predicted for various rtrana and 



/Sciuatcr values in Figure [7]i and b. Significantly, this implies 
that very similar RLFs are also predicted for sources of very 
different ages, so long as the duty cycle (i.e. ton/tott) remains 
the same, and sources are not old enough for inverse Comp- 
ton losses to dominate. This is shown in Figure [T);, where 
we vary the jet on time, while keeping ton/^off constant. The 
duty cycle is the main factor determining the number of ra- 
dio quiet sources, since the fraction of sources old enough to 
drop below the luminosity detection threshold is expected to 
be relatively small. This is illustrated in Figure[7li, where du- 
ration of the active phase is kept constant, while the length 
of the quiescent phase is varied. 

The flatness of the tracks for t < few Myrs for a range of 
power- law exponents in Figure [6] also implies that the RLF 
is also not very sensitive to changes in the density profile 
within the galaxy. This is seen in Figure [Tjs. The "charac- 
teristic" luminosity of a source before it suffers significant 
inverse Compton losses is largely determined by mean jet 
power (Figure [Tlf), core radius (Figure [T)?) and density (Fig- 
ure [Tfi), and axial ratio (Figure O) of the source, in the 
sense that higher values of Qjot, rcore, Pcore and lower value 
of Rt result in higher characteristic luminosities. Finally, 
effects of scatter in jet power are shown in Figure [7];. As 
expected, larger scatter in jet power results in less sources 
at the break (or "characteristic") luminosity, and hence a 
broader distribution of sources across the radio luminosity 
bins. The crucial feature of these plots is that they clearly 
show that changes in most of the parameters only result in 
redistributing the radio loud sources in luminosity. The only 
parameter that can significantly alter the predicted radio 
loud fraction is toff/ton, i.e. the duty cycle. Hence observed 
radio loud fractions place tight constraints on relative dura- 
tions of the radio active and quiescent phases in our sample. 

Our sample is complete to L1.4 = 8 x 10^^ WHz~^, 
rather than the 10^^ WHz~^ plotted in Figure [7] There- 
fore, only sources with luminosities brighter than this value 
can be used to constrain the models. Inspection of Figure [7] 
shows that the discussion above is still applicable, although 
jet and environmental parameters affect the total radio loud 
fraction to a larger degree for this higher L1.4 value. 



5 DURATION OF RADIO ACTIVE AND 
QUIESCENT PHASES 

5.1 Source sizes 

It is clear from Section r4.2l that observed RLFs can be used 
to constrain the duty cycles (i.e. relative lengths of the radio 
active and quiescent phases), however it is difficult to place 
constraints on actual values of tan and toft- Instead, only 
tentative upper or lower limits can be placed on these. The 
former would correspond to a case where the best fit to the 
observed RLF is for a population of sources that do not 
suffer significant inverse Compton losses; and the latter to a 
case when they do. 

A much better constraint on these timescales is ob- 
tained by considering source sizes. For a constant jet power, 
dynamical models of Kaiser & Alexander (|1997|) . Kaiser 
et al. (|1997p and Alexander (|2000p predict source sizes as 
well as radio luminosities as a function of time (see Figure|6)). 
Sources are observed at various stages in their lifetimes, and 
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Figure 7. Dependence of the cumulative radio luminosity function (RLF) on input parameters. Base quantities used arc: Qj^t = 
3 X 1035±0-8 W, ncorc = 0.2 cm-3^ ^corc = 1.0 kpc, /3gaiaxy = 1.0, /3ciu=tor = 1-9, Ttrans = 50 kpc, cocoou axlal ratlo _Rt = 2.0, 
ton = 2 X 10* years and fott = 4 X 10* years. Transition radius rtrans between galaxy and cluster- type atmospheres (o), cluster power-law 
atmosphere exponent /3ciuator (fe), duration of on phase ton (c), duty cycle ton/totf {d) and power-law exponent within the galaxy /3gaiaxy 
(e) are not crucial to the shape of the RLF, with the notable exception of the duty cycle determining the number of radio quiet sources 
(see text for a discussion). The distribution of radio loud sources in luminosity is largely determined by mean jet power Qj^t if) and 
scatter (j) , core radius rcore (fi), core density ricore (h), and the cocoon axial ratio -Rt («)• 
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it is the oldest sources that give an estimate of either the 
duration of the active phase, or the age of the source when 
its radio luminosity falls below our detection threshold. The 
discussion of Section 14.11 suggests the latter case will only 
occur for very weak jets, and even then this is unlikely. 

Contour maps of the 1191 radio loud sources in the 
present sample allow source sizes to be determined. Sizes 
of sources with prominent lobes (FR-IIs) were determined 
manually. Defining source size for low luminosity FR-Is 
is not easy, since the source surface brightness decreases 
smoothly with distance from the nucleus. For the purposes 
of a comparative analysis within the presented sample, how- 
ever, it is sufficient to adopt catalogued NVSS and FIRST 
major and minor axes FWHM. 

Strictly speaking, the radio source model employed in 
this work only describes the evolution of FR-II sources. How- 
ever, two points justify its applicability to the whole radio 
source population in our sample. Figure [T^ shows that there 
are relatively few bona fide FR-Is. The largest sources are 
invariably FR-IIs, and at smaller sizes the sample is domi- 
nated by sources that are only just resolved (and hence are 
those sources for which morphologies cannot be reliably de- 
termined). As all radio sources undergo an initial phase of 
supersonic expansion, we expect the FR-II model to pro- 
vide an accurate description of this sub-population. Once 
the supersonic jets are disrupted, observed distributions of 
source sizes ([Kaiser fc Best 2007|l and theoretical consider- 
ations ([Alexander 2000ll suggest source luminosities quickly 
drop below the detection threshold. We performed the analy- 
sis outlined below on both the full radio source sample, and 
with FR-Is excluded, and obtained identical results. This 
justifies our application of the radio source model to FR-I 
sources, and in what follows we present the results for the 
whole population. 



5.2 Orientation effects 

No spectral information is available for the catalogued radio 
sources. This means various orientation effects must be con- 
sidered when interpreting source properties such as size and 
luminosity. The observed quantities will only correspond to 
actual source parameters for sources that are observed face- 
on. Any other orientation (i.e. beaming) conspires to de- 
crease the apparent source size, and increase its luminosity. 
In other words, the source must be moved right and down 
in the P-D plane from the observed position. 

Consider a source subtending angle cj) to the line-of- 
sight. If the source has size D, its projected size will be 
-Dapp = Dsmcj). The apparent luminosity is related to in- 
trinsic luminosity via ( [Blandford fc Konigl 1979| Lapp = 



L5i 



where ~ 



(1-/3^)''"(1-/3dcos0) 



is the 



Doppler factor for /3d = 
of the source. 

Observations of powerful 
l|Hardcastle et al. 1998|l suggest Usourcc 
adopting q = 0.8 we have 



(1 - 0.8cos(^)" 



/c and a is the spectral index 



radio 



0.8c, 



laxies 
and 




0.4 0.6 

Figure 8. Probability of detecting a source with flux below the 
catalogue threshold, S < Smin- The non-relativistic extended 
component of the flux is not affected by beaming. 



The solid angle subtended by (j) is AQ = 27r(l — cos0). 
In accordance with AGN unification models there should 
be no preferred direction for the beaming, and hence this 
solid angle is directly related to the probability of the ma- 
jor axis of the radio source being directed within of the 
line-of-sight. A convolution of this probability with the ex- 
pression for Lapp in Equation[T]then gives the probability of 
a source with luminosity L (or flux S) below the flux thresh- 
old (3.4 mjy for the present sample) being detected due to 
the beaming. Figure |8] plots this probability as a function of 
intrinsic source flux. 

Most radio sources contain a core component as well as 
more extended emission such as lobes. Since only the rel- 
ativistic core component will be beamed, the level of sam- 
ple contamination by low luminosity sources is sensitive to 
the adopted core-to-total flux ratio, with core-dominated 
sources prone to substantial beaming. In general this ef- 
fect is difficult to quantify. However, the presented sample is 
volume-limited. This means that the objects were selected 
according to their optical properties, which are unaffected 
by orientation. Moreover, the correlation between optical 
and radio luminosities at the low luminosity end of the 
sample (Section 13.11 and Figure [3| suggests there are very 
few sources with uncharacteristically high radio luminosi- 
ties. Thus the radio luminosities of objects in the sample 
are unlikely to be greatly affected by beaming. 

The apparent source sizes will not be affected by the 
beaming by anywhere near the same amount as the lumi- 
nosities. For Ar2/47r = 0.01 Equation [T] gives -Dapp/O ~ 7, 
while I/app/Lcxtondcd = 1/200. For an appreciable solid an- 
gle of Af2/47r = 0.3 this value drops to Dapp/I? ~ 1.4, while 



Lapp/ Lc 



— 1/25. Hence any orientation-related cor- 



Da 



D 



(1) 



rections are likely to affect the derived radio luminosities to a 
much greater extent than source sizes. As a zeroth-order ap- 
proximation, the source will move (almost vertically) down 
in the P-D plane. This introduces larger errors in the derived 
jet powers than source ages. 
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Figure 9. Distribution of jet power as a function of bulge mass. 
These are derived by combining observed source sizes and lumi- 
nosities with the radio source model. Radio source and environ- 
mental parameters of Section 14.11 are used for the modeling. 
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Figure 11. Distribution of jet power with [Olll] line luminosity 
for the radio AGN subsample. The Rawlings & Saunders (1991) 
relation is shown by a dashed line. 
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Figure 10. Fractile distribution of jet power for the top four 
stellar mass bins. Lowest mass bins are not shown due to large 
numbers of unresolved sources in these affecting the statistics 
significantly. Massive galaxies host more powerful jets. 

6 DERIVED SOURCE PROPERTIES 
6.1 Jet powers 

For a given atmosphere and cocoon axial ratio, a source with 
specified radio luminosity and size can only be described 
by a single evolutionary track. Hence, adopting the density 
profiles of Section 14.21 and once again taking J?t = 2, the 
position of a source in the P-D plane uniquely defines its 
age and jet power. Proceeding in this fashion, we derived 
jet powers and ages for every source classified as a radio 
loud AGN. The corresponding jet powers are plotted as a 
function of stellar mass in Figure [S] Although the scatter is 
large, the more massive hosts contain less sources with lower 
(log Qjot < 35.5 W) jet powers. 

This point is illustrated in Figure IIUI where we plot 
the jet power fractile ranges for various stellar mass bins. In 
the bulk of the sample (0.3 to 0.8— fractile ranges) all mass 
bins follow a similar trend, with the more massive hosts 
consistently hosting more powerful radio sources. 

The mean jet power (Table |4ll increases with stellar 



mass as Qjot oc Af*° *=^° ^ Benson et al. (2007) analysed 
a sample of ^ 9000 SDSS galaxies and found that almost 
all stellar mass is located in the spheroid component for 
the stellar masses considered in the present sample. Hence 
Qjet oc Ml^^^-^, or Qjet oc Mbh°-^=^ ° ' using the 0.3 de x 
scatter in the A/bh — Mbuige relation ( [Haring fc Rix 2004p . 
This is significantly less steep than the Qjct oc Afen ex- 
pected from Bondi accretion, suggesting there exists some 
factor limiting the accretion/outflow process. The more mas- 
sive galaxies also show a larger fraction of high Qjc values 
(i.e. a larger scatter at the high-Qjot end; see Table|3|. Since 
it is the number of powerful sources that determines the 
shape of the radio luminosity function at the bright, low 
number counts end, these values for the scatter in jet power 
are adopted when fitting the observed bivariate luminosity 
function below. 

6.2 Independence of optical and radio AGN 
activity 

Best et al. (2005b) considered the relation between emission 
line and radio AGN activity in their sample in two ways. 
First, they looked at the distribution of [QUI] luminosities 
as a function of radio luminosity, and found there to be no 
correlation between the two quantities for their (mostly low 
ii.4GHz) sample (their Fig. 5). Secondly, Best et al. com- 
pared the radio loud AGN fractions as a function of stellar 
mass in emission line and optically inactive AGNs, finding no 
significant difference between the two (their Fig. 7). These 
findings suggest that whether or not a given object is clas- 
sified as an emission line AGN or a normal galaxy has no 
bearing on how likely it is to host a radio AGN. An identical 
analysis of our sample confirmed these results. 

Using radio source models, we extend this analysis to 
investigation of the radio jet power dependence on [OIII] line 
luminosity, as shown in Figure [TT] 

The [OIII] luminosity probes the narrow-line region, and 
is thus indicative of the AGN bolometric luminosity. No cor- 
relation between Qjet and L[oni] is observed in the largely 
low radio luminosity sample presented here. By constrast, 
using spectral ages to estimate jet powers, Rawlings & Saun- 
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Table 4. Jet powers and timescales determined from observed RLFs and source distribution in the P-D plane. Jet powers and median 
source ages are determined from observed source sizes using parameters of Scction l4.2l Quiescent phase duration is found from ratios 
obtained by fitting to the observed luminosity function for each Mi, bin. 



ders (1991) found a tight relation between narrow-line opti- 
cal luminosity and jet power in powerful radio sources; this 
is shown in Figure [Tl]by a dashed line. The fact that no such 
correlation exists for low luminosity radio sources confirms 
the claim of Best et al. that radio and optical AGN activity 
are distinct phenomena in these objects. An intriguing pos- 
sibility is that these correspond to different accretion states. 
In this picture, powerful radio jets would be fuelled when 
the accretion disk is in a radiatively inefficient, low luminos- 
ity state; while optical AGN are detected when the disk is 
radiatively efficient and correspondingly the jet power is low 
dNarayan et al. 1998[rMeier 200ip . 

6.3 Active timescales 

The distribution of sources in the P-D plane can also be 
used to place constraints on radio source ages. Figure 1121 
shows the observed distribution for each M* bin. Also plot- 
ted are P-D tracks for log Qjot = log Qjct ± f7"iogQ using the 
parameters of Section 14.21 and the derived distribution of 
jet powers given in Tabled Curves corresponding to source 
ages of 10* and 3 x 10* years are given for guidance; these 
are the dotted, almost vertical, lines. Finally, limits on the 
maximum detectable size of a source with given luminosity 
and redshift are plotted as straight coloured lines in log-log 
space. Very dim extended sources, such as radio relics, would 
fall below these limits. As discussed in Section [2]2] these are 
largely missed by the pairing process employed in construct- 
ing the present catalogue. However, since in this work the 
emphasis is on the currently active radio sources, this se- 
lection effect does not alter the analysis. More importantly. 
Figure [T2] shows that there are no resolved sources (i.e. those 
sources whose positions are well-defined in the P-D plane) 
close to the NVSS detection threshold (green line). Thus it 
is unlikely that a significant number of large radio sources 
is missed. 

Inspection of Figure [12] shows that although large (and 
hence old) sources are found in all mass bins, the median 
source size is smaller in lower mass bins. Figure [13] quanti- 
fies this statement by using the derived source ages and once 
again plotting the fractile ranges. Unresolved sources pose a 
problem for this approach, since only upper limits on their 
sizes are available. We assume that all unresolved sources 
have very small sizes and are thus younger than the smallest 
resolved sources. While crude, this assumption is sufficiently 
accurate for the highest stellar mass bins, where the unre- 
solved number counts are low, and there are many sources 
larger than the upper unresolved limit. However, no mean- 
ingful information can be extracted for log A^jt/Mq < 11.16. 
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Figure 13. Fractile distribution of radio source ages for the top 
four stellar mass bins. Lowest mass bins are not shown due to large 
numbers of unresolved sources in these affecting the statistics 
significantly. 



Much like jet powers, the derived source ages appear 
to depend strongly on stellar mass, with the most massive 
galaxies hosting older sources. Using the median to„ values 
(Table[4]) we find oc hU^-'^ . 



6.4 Quiescent timescales 

Duration of the quiescent phase fott is found from the bi- 
variate luminosity function. Figure [14] plots the best fits to 
the observed RLF for each stellar mass bin. Parameters of 
Sections 14.21 and 15.11 and jet powers given in Table [3] are 
adopted. The only adjustable parameters are the jet on and 
off times. 

As discussed in Section 14.21 it is difficult to constrain 
the active phase duration from the radio luminosity func- 
tion alone. However, the RLFs are very sensitive to the 
ratio toff /ton- Therefore, combining ton estimates from ob- 
served source sizes with constraints on the duty cycle from 
the luminosity functions yields tott values for each bin. 
These are given in Table [4] where ton is calculated from the 
median values of Table [3] by assuming that a typical radio 
source with a given stellar mass will be active for 2ton, median, 
i.e. it is observed halfway through its evolution. 
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Figure 12. Distribution of radio loud sources in tlic pliysical size - radio luminosity plane for a given range of stellar masses. Crosses 
represent sources whose size is determined by inspection. Open squares represent NVSS, and circles are FIRST sources. Blue points are 
FR-IIs, red FR-Is; and black points are sources that are only just resolved, and hence have uncertain morphologies. Errors on resolved 
source sizes are FWHM. Also plotted are the limits on maximum detectable source size, given the sensitivity of each survey, as a function 
of redshift and luminosity. These are shown as green (NVSS) and purple (FIRST) lines. P-D tracks using parameters of Section |4. 21 with 
jet power corresponding to mean ± cr (Table |4]l are also shown as thick black lines. Finally, we plot the expected tracks corresponding 
to source ages of 10* and 3 X 10*^ years; these are shown by dotted lines, with smaller source sizes corresponding to younger sources. 
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Li.4GHz/WHz-' L,4oHz/WHz-' 
(c) 11.36 < logM*/AfQ < 11.56 (d) 11.16 < logM*/A-/Q < 11.36 



Figure 14. Predicted cumulative radio luminosity functions as a function of stellar mass. The curves represent various jet on times; solid 
= lO'^ yrs, long-dashed = 2.5 X 10^ yrs, short-dashed = 8 X 10^ yrs, dotted = 2 X 10* years, dot-dashed = 6.3 X 10* years. Parameters 
used for fitting are the same as those from which source ages and jet powers were determined. The observed luminosity functions are 
only plotted for L1.4 > 8 X 10^^ WHz~^ as the sample is incomplete below this luminosity. No plot is shown for logM*/M0 < 11.16 
due to large uncertainties in derived jet powers. 



7 DISCUSSION 

7.1 Robustness of timescale estimates 

It is important to consider how reliable the timescales de- 
rived in the preceding section are. The ratio tott/ion ratio 
is largely set by the radio loud fraction in each bin (see 
Section 14. 2|) . and hence uncertainties in this quantity are 
comparable with Poisson errors in the corresponding lumi- 
nosity function. These are clearly negligible in comparison 
with variability in parameters such as the density profile, jet 
power and cocoon axial ratio, all of which are important for 
individual source age determination. 

The dynamical model of Kaiser & Alexander (|1997|l re- 
lates the size of a source expanding in a power-law atmo- 
sphere with exponent /3 to its age t and various cocoon and 
ICM parameters via 



D ocRt: 



t 



(2) 



Here, the cocoon axial ratio Rt is inversely proportional 
to jet opening angle 9. For the majority of our sample we 
do not expect 8 to vary by more than a factor of two, and 
as t oc i?T~*''^ we do not expect the axial ratio to affect our 
source age determinations by more than about a factor of 
two also. 



Neglecting relativistic electron energy loss processes 
(adiabatic and radiative). Kaiser et al. (|1997|) model gives 



o 5±£ 5+p 16-4p-(p + 5)/3 

OC (Poororooro'^) Qjot D 1 = 



(3) 



where p = 2.14 is the power-law exponent of the elec- 
tron energy distribution. Rearranging yields 



(4) 



where x 



and y 



In 



66+6p-3(p-H5))3 

12(5-^) y ^ 12 (5 --3 ) 

deriving individual source jet powers in Section lSTTI we had to 
assume a density profile. However, from Equation[2]it is clear 
that source size does not depend on Qjct or the density pro- 
file individually, but instead the combination q = — — ^ . 
By knowing both the source radio luminosity anif size, and 

( ... , 

using equations [2] and |4] we arrive at t oc g V \^ v / 
Taking (3 — 1.5 gives x = 1.1 and y = 0.6, yielding t oc g""'^. 
Hence even an uncertainty of an order of magnitude in q will 
alter the derived timescale by only a factor of a few. We can 
thus conclude that our tan estimates are likely to be correct 
to within a factor of two or so. As the uncertainty associated 
with the ratio tott/^on is much less than this, duration of the 
quiescent timescale is also expected to be similarly accurate. 
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7.2 Mass dependence 

Best et al. H2005b|) found a strong dependence of the ra- 
dio fond fraction on black hole and host galaxy mass, 
fRL oc Mbh^ '^*° ^ We find a similar result for our volume- 
limited sample (Figurel!]), with Jrl oc Af*^'^*"'^. For a fixed 
spheroid-to-total light ratio (Section 16. ip . using the Mbh — 
relatio n paring fc Rix 2004t gives Jrl oc Mun'^^^-^. 
Best et al. H2005b|) point out that this relation is similar to 
the Mcooi oc Mbh^'^ gas cooling rate dependence on black 
hole mass derived from the X-ray-optical correlation in lumi- 
nous elliptical gala^xies, and suggest that AGN radio activity 
might be fuelled by the cooling of hot gas within the host 
galaxy. 

The derived AGN "on" timescales in Table |4] clearly in- 
crease with stellar mass. If the AGN is fuelled by cool gas 
accretion, and the duration of the active phase is limited 
by fuel availability, one would expect for constant radia- 
tive and accretion efficiency fon oc ^^s°L^ which on using 

^jct 

the scaling relations of Section l6.ll gives ton oc j\^^0-9±0'5^ 
This is consistent with ton oc M*^'^ derived in Section [6.31 
and suggests that the duration of the active phase may in- 
deed be determined by the availability of fuel to "feed" the 
jets. The quiescent phase duration (Table [4]) also shows a 
strong dependence on Mi^, with sources hosted by the most 
massive galaxies spending a shorter amount of time in the 
off phase, consistent with higher cooling rates in those ob- 
jects. In the absence of a heating source one would expect 
tott oc Mcooi^^ oc M*"^'^*"'^. However, in practice as the 
more massive galaxies host the more powerful radio sources, 
they will be more affected by the heating and cold gas up- 
lifting than their lower mass counterparts, and the observed 
relation is not as steep. 



8 CONCLUSIONS 

We constructed a fiux- and volume-limited (0.03 < z < 0.1) 
sample of radio sources with optical identifications by cross 
correlating the SDSS optical survey with 1.4 GHz NVSS and 
FIRST surveys. Source sizes and luminosities together with 
radio source models allowed us to derive jet powers and ages 
of individual sources, and hence the jet "on" time as a func- 
tion of stellar mass. The bivariate luminosity function was 
then used to constrain the time a typical radio source spends 
in an inactive state. Radio and emission line AGN activity 
are found to be independent phenomena. We also find that 
both the radio source lifetime and duration of the quiescent 
phase have a strong mass dependence, with massive hosts 
harbouring longer-lived sources that are triggered more fre- 
quently. Gas cooling rate shows a similar mass dependence, 
suggesting that fuel depletion is the reason the jets switch 
off. 
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